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EXECUTIVE  SUMMARY 


Introduction 

The  Defense  Communications  Agency  (DCA)  has  been  evaluating  the  effectiveness  of  horizontal 
electric  polarized  waves  for  high-altitude  air-to-air  communication.  Horizontal  polarization  is  the  best 
choice  because  a  high-speed  high-altitude  jet  aircraft  can  trail  a  horizontal  long  wire  antenna,  which  pre¬ 
ferentially  excites  transverse-electric  (TE)  radio  waves,  and  because  the  ambient  noise  caused  by  the 
vertically  polarized  lightning  strokes,  which  is  most  of  the  ambient  noise  in  VLF  radio  systems,  is 
believed  to  be  least  at  horizontal  polarization.  A  computer-program  model  has  been  developed  for 
predicting  the  horizontally  and  vertically  polarized  atmospheric  noise  at  any  location,  season,  time  of 
day,  and  altitude.  This  noise  model  will  play  a  role  in  the  prediction  of  VLF  system  performance. 

The  Program  For  The  Atmospheric-Noise  Model 

The  program  for  the  atmospheric-noise  model  combines  two  previously  written  programs  called 
COMPWR  and  NOISLAN  with  a  new  program  called  HORNS.  The  COMPWR  and  NOISLAN  pro¬ 
grams,  produced  oy  the  Westinghouse  Georesearch  Laboratory  (WGL)  under  contract  to  NRL,  predict 
the  vertically  polarized  electric  atmospheric  noise  in  the  frequency  band  from  10  to  30  kHz  in  the 
quasi-transverse-magnetic  (quasi-TM)  waveguide  modes  at  the  surface  of  the  earth  at  any  location,  time 
of  day,  or  season.  The  HORNS  program  calculates  the  horizontally  and  vertically  polarized  electric  and 
magnetic  fields  at  any  altitude  in  the  earth-ionosphere  waveguide,  when  the  vertical  electric  fields  in  the 
propagating  modes  at  the  earth’s  surface  are  known.  The  predictions  of  the  COMPWR  and  NOISLAN 
programs  compare  well  with  extensive  measurements  of  worldwide  atmospheric  noise.  The  predictions 
of  the  algorithms  of  the  HORNS  program  for  the  high-altitude  fields  have  been  compared  with  experi¬ 
mental  data  to  some  extent,  and  they  should  be  reasonably  accurate,  since  they  are  based  on  full-wave 
solutions  of  Maxwell’s  equations.  The  COMPWR  program,  which  calculates  the  vertical  effective  radi¬ 
ated  power  (VERP)  for  each  noise  source,  has  two  available  options.  COMPWR  option  I  is  the  original 
WGL  model.  COMPWR  option  II  contains  different  noise-source  powers  and  diurnal  power-modifier 
functions  that  vary  the  VERP  of  each  atmospheric  noise  source  over  a  24-hour  period.  The  NOISLAN 
program,  which  calculates  propagation  effects  ant*  sums  the  field  strength  contributions  from  all  the 
sources,  also  has  two  available  options.  Option  I  is  the  original  WGL  model.  Option  II  contains 
isotropic-waveguide-model  attenuation  rates  and  excitation  factors  as  computed  by  a  program  called 
WAVEGUID  using  exponential  ionosphere  profiles  with  ft  =  0.3  km-1  and  h  =  70  km  for  day  and  with 
li  =  0.5  km” 1  and  li  =  90  km  for  night.  The  HORNS  program  calculates  components  of  the  high- 
altitude  electromagnetic  noise  field  propagated  by  three  quasi-TM  modes  in  the  earth-ionosphere 
waveguide.  The  noise  propagated  in  quasi-TE  modes  is  not  evaluated. 

Comparison  With  Data 

High-altitude  noise  data  available  Tor  analyzing  the  HORNS  model  are  limited.  The  Mitre  Cor¬ 
poration  has  made  a  few  airborne  noise  measurements  suitable  for  comparison  with  tl  e  atmospheric 
noise  predicted  by  the  HORNS  computer  program.  Reasonably  good  agreement  exists  between  these 
measured  values  and  the  predicted  values  of  both  the  vertically  and  horizontally  polarized  noise,  with 
almost  all  the  data  being  within  1  standard  deviation  of  the  predicted  values,  if  the  standard  deviation  is 
calculated  from  the  values  of  given  in  CCIR  Report  322  111.  This  result  is  encouraging,  but  more 
data  are  needed  to  fully  validate  the  model. 


Additional  Development  of  the  Model 


Several  possibilities  exist  for  further  developing  the  model.  However,  the  most  useful  additional 
development  at  present  would  be  to  extend  the  model’s  frequency  range,  now  10  to  30  kHz.  to  60  kllz, 
since  this  extended  frequency  range  is  of  great  interest  for  the  DC  A  MliliCN  (Minimum  Essential 
F.mergency  Communications  Network),  and  of  operational  concern  to  the  Navy  and  the  Air  Force. 
Such  an  extension  would  first  require  that  the  COMPWR  and  NOISLAN  programs  be  modified  to  accu¬ 
rately  predict  vertically  polarized  noise  at  the  earth's  surface  up  to  60  kllz.  The  validity  of  the  predic¬ 
tions  Irom  30  to  60  kHz  could  be  checked  against  available  measured  vertical  noise  data,  just  as  was 
done  for  10  to  30  kllz.  Then,  the  HORNS  program  could  easily  be  extended  to  calculate  high-altitude 
noise.  Another  improvement  in  the  atmospheric-noise  prediction  model  would  be  to  include  the  noise 
contributed  by  the  quasi-TF.  mode. 


COMPUTER-PROGRAM  MODEL  FOR  PREDICTING 
HORIZONTALLY  AND  VERTICALLY  POLARIZED 
VLF  ATMOSPHERIC  RADIO  NOISE  AT  ELEVATED  RECEIVERS 


INTRODUCTION 

The  communication  effectiveness  of  any  radio  circuit  or  link  strongly  depends  on  the  signal-to- 
noise  ratio  in  the  receiver.  At  very  low  frequencies  the  propagated  radio  noise  generated  by  lightning 
discharges,  referred  to  as  atmospheric  radio  noise,  is  frequently  the  dominant  system  noise. 

Many  workers  have  studied  atmospheric  radio  noise  (1-7).  In  particular,  Ref.  3  decribes  a  com¬ 
puter  model  which  predicts  the  atmospheric  radio  noise  generated  by  thunderstorms  in  all  regions  of 
the  world  and  then  calculates  the  propagation  of  the  noise  to  a  receiver  at  any  point  on  the  earth’s  sur¬ 
face.  By  performing  this  computation  over  a  grid  of  receiver  locations,  a  map  of  the  intensity  of  atmos¬ 
pheric  radio  noise  can  be  generated  (8)  similar  to  those  given  in  Ref.  1.  However,  as  with  most  work  in 
this  field,  only  the  vertical  electric  component  of  the  noise  is  calculated. 

In  recent  years  interest  has  grown  in  the  use  of  horizontal  electric  polarized  waves  for  high- 
altitude  air-to-air  communication.  This  polarization  should  be  advantageous,  because  it  is  operationally 
easier  for  a  high-speed  high-altitude  jet  aircraft  to  trail  a  horizontal  long  wire  antenna,  which  preferen¬ 
tially  excites  transverse  electic  radio  waves,  and  because  the  vertically  polarized  lightning  strokes  are 
thought  to  be  less  effective  in  generating  horizontally  polarized  noises.  Therefore,  a  requirement  exists 
to  be  able  to  accurately  predict  these  noise  levels. 

The  theoretical  and  experimental  effects  associated  with  horizontally  polarized  VLF-wave 
transmission  and  reception  have  been  treated  by  several  authors  [9-13],  Also,  the  use  of  airborne,  hor¬ 
izontally  polarized  transmitting  and  receiving  antennas  for  communications  at  VLF  has  been  discussed 
several  times  [14-17],  The  WAVEGUID  computer  program  [11]  contains  the  algorithms  necessary  to 
compute  all  the  electric  and  magnetic  components  cf  a  propagating  electromagnetic  wave  in  a  given 
waveguide  mode  at  any  height  above  the  ground,  if  the  vertical  electric  component  of  the  wave  at  the 
ground  is  known.  The  COMPWR  and  NOISLAN  programs  produced  by  the  Westinghouse 
Georesearch  Laboratory  (WGL)  [3]  calculate  the  amplitudes  of  the  vertical  electric  VLF  noise  in  the 
three  dominent  quasi-transverse-magnetic  (quasi-TM)  waveguide  modes  from  each  equivalent  atmos¬ 
pheric  noise  source.  By  combining  these  programs,  NRL  has  now  developed  a  working  computer- 
program  model  for  predicting  the  horizontally  and  vertically  polarized  atmospheric  noise  at  any  location, 
season,  time  of  day,  and  altitude.  This  noise  program  should  assist  in  the  prediction  of  system  perfor¬ 
mance. 

THE  PROGRAMS  THAT  FORM  THE  NOISE-PREDICTING  MODEL 

References  3  and  11  give  complete  descriptions  of  the  COMPWR,  NOISLAN,  and  WAVEGUID 
programs.  The  following  subsections  briefly  describe  the  features  of  each  program  essential  to  the  new 
program. 
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Description  of  COMPWR  and  NOISLAN 

The  atmospheric-noise  model  is  built  on  thunderstorm-day  data,  from  which  vertical  electrical 
noise  intensity  is  derived,  The  thunderstorm-day  is  a  weather  statistic  which  gives  over  a  30-day  inter¬ 
val  the  number  of  days  during  which  at  least  one  thunderstorm  occurred  in  the  vicinity  of  the  observa¬ 
tion  site.  The  data  base  used  by  the  program  represents  observations  collected  over  several  decades  and 
at  thousands  of  locations.  In  addition  the  model  uses  data  on  the  diurnal  variation  and  variability  of 
thunderstorm  occurrences.  Algorithms  convert  thunderstorm  days  into  the  number  and  type  of  lightn¬ 
ing  discharges  per  unit  area  on  the  earth’s  surface.  Also,  the  current  moments,  waveforms,  and  power 
spectral  densities  of  the  pulses  are  calculated.  The  resultant  source  powers  are  combined  into 
equivalent  noise  "transmitters"  representing  5°  by  5°  areas  of  the  earth’s  surface.  These  5°  by  5°  areas 
are  further  combined  into  15°  by  15°  areas.  The  program  decides,  based  on  the  proximity  and  power  of 
the  equivalent  noise  transmitters,  whether  to  use  5°  by  5°  or  15°  by  15°  areas  in  the  final  computations. 


After  COMPWR  establishes  the  location,  power,  and  standard  deviation  of  the  power  for  each 
equivalent  noise  transmitter,  the  NOISLAN  program  then  computes  the  resulant  vertical  electric  field  at 
the  specified  receiver  location  using  the  propagation  equation 
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is  the  vertical  electric  component  of  the  atmospheric  noise  contribution  of  the  Ah  noise 
source  in  the  mh  waveguide  mode  at  the  surface  of  the  waveguide,  in  volts  per  meter, 
is  a  normalizing  constant  equal  to  1.58  x  105  , 

is  the  vertical  effective  radiated  power  (VERP)  of  the  Ah  equivalent  noise  transmitter 
in  watts, 

is  the  excitation  factor  of  the  mh  waveguide  mode  at  the  Ah  transmitter  locations, 
is  the  excitation  factor  of  the  nth  waveguide  mode  at  the  receiver  location, 
is  the  frequency  in  hertz, 

are  the  effective  ionospheric  heights  in  meters  at  the  Ah  transmitter  and  at  the  receiver 
respectively, 

is  the  earth’s  radius  in  meters, 

is  the  path  length  from  the  Ah  transmitter  to  the  receiver  in  meters, 


and 

a,'f)  is  the  effective  attenuation  rate  in  decibels  per  megameter  (1000  km)  for  the 

mh  mode  from  the  Ah  transmitter. 


The  effective  attenuation  rate  a,*/’  is  computed  from  the  combination  of  several  attenuation  terms  as 
follows: 
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In  Eq.  (2),  '  is  a  directional  factor  and  is  given  by 
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where  is  the  direction  of  propagation  of  the  wave  from  the  Ah  transmitter  with  respect  to  geo¬ 
graphic  north  and  k  is  given  by 
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in  which 

0  is  the  colatitude  of  the  receiver  in  degrees, 

/  is  the  frequency  in  hertz, 

Km is  a  constant  that  for  west-to-east  propagation  equals  0.68  during  the  day  und  0.70  at  night 
and  for  east-to-west  propagation  equals  1.10  during  the  day  and  2.00  at  night, 

and 

(3  is  a  constant  that  equals  0.04  during  the  day  and  0.03  at  night. 

Also  in  Eq.  (2),  is  the  attenuation  rate  (day  or  night)  of  the  mh-order  mode  over  seawater 
for  the  isotropic  case,  A«  jjn  is  the  differential  ionospheric  attenuation  rate  in  effect  for  the  incremental 
distance  Ad(  along  the  path,  Aa  is  the  differential  ground  attenuation  rate  due  to  ground  conductivi¬ 
ties  other  than  seawater,  and  A <4  is  the  distance  segment  of  the  propagation  path  within  a  given  con¬ 
ductivity  area.  The  differential  ionospheric  attenuation  term  is  obtained  from 

A«  u„  **  fl2A  cos +  £(sin  lx /I)  sin^,],  (5) 

where 

0 1  is  the  colatitude  of  the  Ah  segment  of  the  propagation  path  in  degrees, 

Xi  is  the  solar  zenith  angle  of  the  Ah  segment  of  the  propagation  path, 

A  equals  2.2  x  10~5  dB/Mm-Hz, 

and 

B  equals  2.8  x  10-5  dB/Mm-Hz. 

The  values  in  Eq.  (2)  for  A  a  £/{.  ,,  are  the  differences  between  the  attenuation  rate  of  seawater  for  the 
isotropic  case  and  the  attenuation  rates  for  the  lower  ground  conductivities  as  stated  in  the  equa¬ 
tion 

-  "cl. L 

where  values  of a^j\„  andaa„  are  those  given  by  Wait  and  Spies  [181  for  h  =»  70  km  and/3  -  0.3  km-1 
for  day  and  h  -  90  km  and  /3  -  0.5  km'1  for  night  with  11  -  0  (with  their  parameter  not  being  the 
same  as  the  parameter  B  in  Eq.  (4)).  However,  Ref.  18  is  incomplete,  because  no  theoretical  data  were 
obtained  for  some  of  the  lower  values  of  ground  conductivity  for  the  secor.d-order  and  third-order 
modes.  When  the  theoretical  data  were  not  available,  the  authors  of  the  NOISLAN  program  appear  to 
have  extrapolated  data  according  to  an  unspecified  procedure.  The  values  of  Ad*  are  computed  from  a 
digitized  version  [3,  pp.  4-191  of  a  map  of  VLF  ground  conductivity  developed  for  NRL  by  Morgan 
119], 


The  field-strength  contributions  from  each  transmitter  and  each  mode  are  then  summed  in  an 
RMS  manner,  thus  giving  the  vertical  electric  field  strength  £;(0)  of  the  noise  at  the  ground  for  a 
specific  hour  of  the  day,  month,  location,  and  frequency; 
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NOISLAN  also  computes  the  standard  deviation  of  the  noise  <r„  and  the  voltage  deviation  f',,.  How 
ever,  a  discussion  of  these  two  parameters  is  not  pertinent  to  this  report. 


Additional  Options  for  COMPWR  and  NOISLAN 

Two  options  have  been  added  to  the  WGL  model.  The  first  involves  a  change  in  the  part  of  the 
WGL  model  contained  in  program  COMPWR,  The  original  version  of  COMPWR  has  been  designated 
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as  option  1  of  the  COMPWR  program,  and  the  modilkation  to  be  described  has  been  designated  as 
option  11.  COMPWR  computes  the  average  vertical  effective  radiated  power  (VERP)  and  the  standard 
deviation  of, the  VERP  for  each  equivalent  noise  transmitter  for  each  month  of  the  year.  COMPWR 
then  multiplies  the  average  VERP  by  a  diurnal  modifier  to  give  the  VERP  for  each  hour  of  the  day. 
Option  II  involves  changes  to  the  diurnal  modifiers  as  well  as  inclusion  of  an  additional  set  of  modifiers 
to  adjust  the  average  VERP.  These  changes  were  developed  empirically  120,21], 

Another  option  has  been  added  in  the  form  of  a  modification  to  the  NOISLAN  program.  The  ori* 
ginal  version  of  NOISLAN  has  been  designated  us  option  I,  and  the  new  modification  has  been  desig¬ 
nated  as  option  II.  The  modification  involves  the  terms  A,,,  and  A },'}  from  Eq.  (1)  and  and 
(V1, „  from  Eq.  (2).  In  option  I  of  the  NOISLAN  program  the  values  for  these  terms  are  those  given 
by  Watt  [221,  who  in  turn  derived  the  values  from  the  work  of  Wait  and  Spies  [181.  However,  the 
computer  progrum  used  to  generate  the  data  of  Ref.  18  (in  many  cuses)  could  find  no  solution  for  the 
second-order  and  third-order  modes,  especially  for  the  lower  values  of  ground  conductivity,  Further, 
the  accuracy  of  Wait’s  second-order-mode  and  third-order-mode  roots  has  been  questioned  by  Pappert 
et  al.  [23] .  Because  of  these  problems,  it  was  decided  to  recompute  the  appropriate  mode  parameters 
using  the  WAVEGUID  program.  These  parameters  were  then  used  to  compute  the  aforementioned 
terms  from  Eqs.  (1)  and  (2),  The  values  of  the  propagation  parameters  computed  from  the 
WAVEGUID  program  agreed  with  the  values  used  in  the  NOISLAN  program  for  the  first-order  mode, 
but  disagreed  (as  noted  by  Pappert  [231)  for  the  second-order  and  third-order  modes.  Since  the 
second-order  and  third-order  modes  often  carry  more  horizontally  polarized  energy  than  does  the  first- 
order  mode,  accurate  computation  of  the  higher  order  modes  is  important  in  the  prediction  of  high- 
altitude  atmospheric  radio  noise.  The  WAVEGUID  mode-parameter  solutions  are  believed  tc  be  more 
accurate  than  those  of  Refs.  18  and  22.  Option  II,  therefore,  has  been  included  in  NOISLAN  to  use 
the  newly  calculated  parameters.  Also  there  has  been  interest  in  extending  the  VLF  atmospheric  noise 
model  upward  in  frequency  to  cover  the  lower  end  of  the  LF  spectrum  (30  to  60  kHz).  This  extension 
would  not  be  possible  using  the  results  of  Refs.  18  and  22  but  would  be  possible  by  use  of  the 
WAVEGUID  program. 

Description  of  WAVEGUID 

The  WAVEGUID  computer  program  was  originally  developed  at  the  Naval  Electronics  Laboratory 
Center  to  predict  the  vertical  electric  field  received  at  a  point  on  the  earth's  surface  produced  by  a  verti¬ 
cal  electric  transmitting  antenna  at  another  point  on  the  earth’s  surface.  It  was  later  extended  to  give 
the  crosspolarized  (horizontal)  components  of  the  field  at  any  height  in  the  earth-ionosphere  waveguide 
and  even  in  the  ionosphere  itself  [1  1,  24,  25].  The  basic  formulation  of  the  equations  for  the  fields  are 
given  in  Ref.  13.  The  net  resultant  vertical  electric  field  at  the  earth’s  surface  is  represented  as  a  sum¬ 
mation  of  waveguide-mode  fields: 

£•  (0)  -  1  Ids  -v/  ~  dJ?r Y  £[sinVJtf„  exp(/ir/4  +  ik0ds\nQ „)],  (7) 

V<A  h  V  sin(a/a)  „ 

where 
Ids 


a 
d 
h 

A'o 

\ 

and 

Mo  and  € o 

4 


is  the  eigenangle  of  the  mh  waveguide  mode, 

is  the  dipole  moment  (ampere-meter)  of  the  vertical  antenna  located  on  the  earth’s  sur¬ 
face, 

is  the  excitation  factor  of  the  mh  waveguide  mode, 
is  the  radius  of  the  earth, 

is  the  great-circle  distance  between  the  transmitter  and  the  receiver, 
is  the  reference  height  of  the  ionosphere  (used  in  the  definition  of  \„), 

=  co/e  is  the  free-space  propagation  constant  of  the  wave, 
is  the  wavelength  of  the  wave, 

are  the  magnetic  permeability  and  dielectric  constant  of  free  space. 
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Each  waveguide  mode  is  charueieri/.ed  by  its  own  value  of  attenuation  rate  and  phase  velocity 
i '//,  which  are  related  to  the  eigenangle  of  the  «th  waveguide  mode  0 „  according  to 

0.02895  o>  lm (sin  «„)  (8) 

and 


''  Re(sin«M)  ' 

where  « „  is  measured  in  decibels  per  1000  km  (dB/Mm)  and  c  is  the  speed  of  light. 

The  values  of  the  waveguide-mode  eigenangles  depend  on  the  reflection  coefficients  of  the  ground 
and  the  ionosphere.  Because  of  its  anisotropy,  the  ionosphere  hus  four  reflection  coefficients  M  /? n (w  „) , 
l  R[  (0„),  ii  /?,  (W  ,  and  i  Ru(<>  „)  for  a  given  angle  of  incidence  0  The  ground  is  assumed  to  be  isotro¬ 
pic  and  to  have  reflection  coefficients  M /fM (rt H)  and  i  Rl^h)<  which  are  calculable  from  the  eigenangle, 
ground  conductivity,  and  dielectric  constant.  The  WAVEGUID  program  uses  a  procedure  for  calculat¬ 
ing  the  ground  and  ionospheric  reflection  coefficients  referenced  to  any  height  z  in  the  waveguide. 
References  26,  27,  and  28  discuss  the  reflection  coefficients  further.  The  eigenangle  0 „  for  a  given 
waveguide  mode  is  obtained  by  satisfying  the  mode  equation 
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This  matrix  equation  is  simply  rewritten 

t  (o„)  —  (|i flu  nfl||- 1)  i  Ri  ~ I)  ~  i  Rn  i  Rin^i  ii/?u  *  0,  (11) 

where  the  dependence  of  the  reflection  coefficient  on  0„  has  been  suppressed  to  simplify  the  notation. 
The  new  variable  /■’(#„)  is  defined  also.  Each  value  of  0„  which  satisfies  Eq.  (11)  is  the  eigenangle  of 
the  mh  waveguide  mode.  The  excitation  factor  \„  for  the  nth  waveguide  mode  is  obtained  from 

.  .  kh  .  (1  +  ii^ii)2  (1  - 1  R\  i  R\ )  .... 


.  kh  .  n 
-i  -y  sin 


where  the  variable  F(f>)  is  as  defined  in  Eq.  (II).  In  addition  to  the  vertical  electric  field  E.  for  a  given 
waveguide  mode  having  eigenangle  „,  in  general  there  exist  five  other  nonzero  field  components  Ex> 
£,,,  HXs  //,. >  and  H.  within  the  waveguide  for  each  mode. 

These  extra  fields  arise  because  the  ionosphere  is  anisotropic  on  account  of  the  earth’s  magnetic 
field.  Propagation  through  and  reflection  from  such  an  anisotropic  medium  rotates  the  plane  of  polari¬ 
zation  of  the  incident  wave  to  generate  crosspelarized  reflected  fields.  Alternatively  one  could  consider 
the  linearly  polarized  wave  transmitted  from  the  vertical  electric  dipole  antenna  to  be  a  linear  superposi¬ 
tion  of  right-handed  and  left-handed  circularly  polarized  waves,  each  reflected  with  a  different  ampli¬ 
tude  at  the  ionosphere.  The  ionospheric  reflection  causes  a  mixture  of  wave  polarizations  at  a  reception 
point  on  or  above  the  earth’s  surface.  The  fields  at  height  z  are  related  to  the  vertical  electric  fields  at 

the  earth’s  surface  £.(0)  in  a  given  waveguide  mode  by  the  following  functions  of  reflection  coefficients 

and  Hankel  functions: 

EXzMEm  -  71,  (z).  (13) 

Ex(z)/  /:.(())  =  ,c(e)/.S.  (14) 

(0) - JttU+JVAU) -  „5, 
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In  these  expressions 


D 


and 


k 

S  and  C 

h\(q)  and  h2(q) 


is  the  radius  of  the  earth, 

is  the  height  of  the  receiver  above  ground, 

is  the  index  of  refraction  of  the  ground  surface, 

is  the  reference  altitude  in  the  waveguide  for  evaluating  the  mode  parameters  and  the 
reflection  coefficients, 

is  the  propagation  constant  of  the  waves  in  free  space  (  =  t»/c), 
are  the  sine  and  cosine  of  the  eigenanide  0„, 
are  the  modified  Hankel  functions  of  order  1/3, 


T) 


is  the  impedance  of  free  space. 


Using  these  relations,  one  can  obtain  a  desired  field  of  a  given  mode  at  any  altitude  z  if  one 
knows  the  vertical  field  in  that  mode  at  the  earth’s  surface.  The  NOISLAN  prediction  program  calcu¬ 
lates  the  vertical  noise  fields  in  the  three  strongest  quasi-TM  waveguide  modes  in  making  its  noise 
prediction.  So,  the  values  of  £.(0)  for  the  three  most  important  modes  are  available,  and  one  can 
obtain  the  crosspolarized  fields  at  any  height  z  by  applying  Eqs.  (13)  through  (18)  to  the  NOISLAN 
£r(0)  values. 


HORNS  Program  and  Approximations 

The  HORNS  portion  of  the  atmospheric-noise  model  predicts  high-altitude  atmospheric  noise  at 
any  polarization  by  using  the  functions  of  the  WAVEGUID  program,  as  summarized  in  Eqs  (13) 
through  (18),  and  the  values  of  £2(0)  computed  by  the  NOISLAN  program.  Equations  (13)  through 
(18)  may  be  rewritten  in  a  notation  more  appropriate  to  the  summation  of  fields  from  many  sources: 


and 


o*)  =  fiftiz)  E“H 0), 

(13') 

o*)  -  gHliHz)E^(0)/Sn{l) 

(14') 

p</>  f ~  \  -iiOl  +i^,)/i(.«(z)£n,',(0) 
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(17') 

««/»!-._  .  dfW(z) 
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(18') 

In  these  expressions  the  field  symbols  such  as  E™(z)  refer  to  the  field  strength  produced  by  the  /th 
noise  source  in  the  mh  waveguide  mode.  The  symbols  for  the  height  gain  functions,  such  as  /„<',>  (r), 
reter  to  these  functions  in  the  «th  waveguide  mode  and  from  the  /'th  transmitter.  Likewise  the  ground 

and  ionospheric  reflection  coefficients  refer  to  particular  modes  and  to  waves  from  a  particular  noise 
source. 
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After  the  individual  contributions  to  the  components  of  the  noise  field  have  been  calculated,  the 
vertically  polarized  electric  and  magnetic  noise  fields  may  be  obtained  by  analogy  to  Eq.  (6): 

A  3  V/; 

I  I  I E"(z)V  (33) 

/-I  «-i 

and 

A  3 

I  I  ■  (34) 

;-i  /i- 1 

A  random-phase  approximation  has  been  made  in  deriving  Eqs.  (33)  and  (34).  In  other  words,  no 
correlations  are  assumed  among  the  elements  E^'J  in  the  sum  of  Eq.  (33)  or  among  the  elements  H ^ 
of  Eq.  (34).  For  different  values  of  /  (different  sources)  the  random-phase  approximation  is  undoubt¬ 
edly  correct,  since  there  is  no  correlation  in  the  arrival  time  of  energy  from  different  sources.  How¬ 
ever,  for  different  values  of  n  (different  modes  from  the  same  source)  the  random-phase  approximation 
is  less  accurate,  since  some  correlation  exists  among  the  arrival  times  of  different  modes  from  the  same 
source.  However,  since  the  noise  "transmitter"  locations  are  only  the  noise  power  centroids  for  noise 
sources  in  their  respective  areas  and  not  the  noise  sources  themselves,  the  actual  propagation  paths 
from  each  lightning  discharge  to  the  receiver  are  not  known.  Thus  one  cannot  accurately  calculate  the 
phase  angles  of  the  various  modes,  leaving  the  random-phase  approximation  as  the  best  alternative. 

To  obtain  equations  for  the  horizontally  polarized  noise  fields,  one  must  first  select  a  horizontal 
axis  direction  .v'  along  which  the  noise  is  to  be  specified.  The  previously  used  coordinates  x,  y,  and  z 
arc  defined  in  relation  to  the  direction  of  propagation  of  a  wave:  z  is  the  vertical  direction,  x  is  horizon¬ 
tal  in  the  direction  of  propagation  of  the  wave,  and  .v  is  horizontal  perpendicular  to  the  direction  of 
wave  propagation.  Since  the  equivalent  sources  are  at  different  bearings  from  the  receiving  location, 
one  cannot  simply  add  the  field  to  that  from  another  source  such  as  H,['x+  11  (z),  because  they 

will  not  be  parallel.  Instead  one  must  find  the  contribution  of  the  various  fields  to  a  resultant  com¬ 
ponent  fixed  in  space  using  the  normal  processes  of  vector  resolution  and  vector  addition.  One  then 
obtains  the  total  horizontal  fields  in  the  fixed  x'  directions  as  follows: 

A3  .1  " 

Ex(z)  =  £  X  I  £«.*•(*)  cos  (fv'v  4-  E^iz)  sindi'v’l2  (35) 

/  -  i  a  -  i 

and 

A  3 

//,.(*>-  jr  £  |  H&'iz)  cos  ffxy  +  H,[\Hz)  sintf|T  I2  .  (36) 

The  resultant  noise  predictions  Ex(z)  and  Hx(z)  are  real  numbers,  and  the  waveguide-mode  field 
amplitudes  E,[x,  f,,0,.1  .  H,[x,  and  /Y,,0,’  are  complex  phasors.  Hence  one  must  use  the  absolute  magni¬ 
tudes  of  these  contributions  in  Eqs.  (35)  and  (36). 

To  determine  Hj'Hz)  from  E,{,‘H 0)  according  to  Eq.  06'),  one  needs  knowledge  of  the 
eigenangle  of  the  zith  waveguide  mode  that  carries  the  wave  from  the  rth  noise  transmitter,  and  of  the 
ionospheric  and  ground  reflection  coefficients  for  this  propagating  mode.  How  can  one  obtain  H  //,? 
One  could  determine  0  hn  from  the  WAVEGUID  program  using  an  appropriate  anisotropic  ionosphere 
model  by  iteratively  solving  Eq.  (11),  as  is  normally  done  with  this  program.  However,  this  procedure 
is  time  consuming  when  applied  to  noise  predictions,  because  it  ordinarily  requires  several  iterations  to 
settle  on  a  satisfactory  value  oH*,,  .  This  undertaking  could  become  expensive  and  tedious,  since  there 
may  be  as  many  as  70  significant  thunderstorm  atmospheric  noise  sources,  each  having  a  different  bear- 
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ing  from  the  receiver  location.  Since  the  angle  of  arrival  of  energy  from  each  source  would  be 
diflerent,  a  new  calculation  off),,  would  be  required  for  each  case.  An  alternative  method  has  been 
used  in  the  present  version  of  the  program. 

The  approximation  for  0n  in  the  present  program  is  to  substitute  the  eigenangle  0n  which  has 
already  been  calculated  for  the  quasi-TM  mode  with  an  isotropic  ionosphere.  This  procedure  makes  the 
assumption  that  the  field  expressions  in  Eqs.  (13')  thorugh  (18')  are  not  sensitive  to  the  value  b„  . 
This  assumption  was  tested  and  partially  verified  for  a  sample  of  noise  predictions.  Certainly  this 
method  of  selecting  is  arbitrary,  and  it  might  be  improved  in  subsequent  versions  of  the  program. 
Once  is  selected,  it  is  used  for  the  evaluation  of  all  the  factors  that  enter  into  the  computation  of  the 
elevated  fields.  The  geomagnetic  field  strength  and  direction  at  the  receiver  location  is  calculated  using 
a  dipole  model  of  the  earth’s  field.  The  Ah  azimuth  angle  required  to  evaluate  the  ionospheric 
reflection  coefficients  is  evaluated  using  this  geomagnetic  model  and  a  knowledge  of  the  direction  to  the 
Ah  equivalent  noise  source. 

The  atmospheric-noise  model  calculates  the  high-altitude  electromagnetic  fields  propagated  in 
three  quasi-TM  modes  in  the  earth-ionosphere  waveguide.  The  fields  generated  by  the  vertical  com¬ 
ponent  of  the  lightning  stroke  bui  propagated  in  the  quasi-transverse-electric  (quasi-TE)  modes  are 
ignored.  This  deficiency  of  the  prediction  model  may  be  remedied  in  future  versions.  It  seems  reason¬ 
able  that  the  contributions  by  the  quasi-TE  mode  to  the  horizontal  electric  and  vertical  magnetic  field 
should  be  comparable  to  the  contributions  by  the  quasi-TM  mode;  hence  the  observed  noise  might 
exceed  the  predicted  values  by  3  dB. 

The  atmospheric-noise  model  also  ignores  the  contributions  of  the  horizontal  component  of  the 
lightning  discharge,  because  it  has  been  assumed  that  the  horizontal  discharges  are  much  weaker  than 
the  vertical  discharges.  Recent  detailed  analyses  of  lightning  discharge  paths  129]  suggest  that  a 
substantial  elevated  horizontal  component  of  the  discharge  channel  existed  in  several  thoroughly 
analyzed  cloud-to-ground  discharges.  This  strong  horizontal  component  would  be  inconsistent  with  our 
present  noise  model  and  suggests  an  important  topic  for  further  investigation.  A  detailed  statistical 
study  of  lightning-channel  structures  and  of  received  noise  polarizations  would  be  useful  in  assessing 
this  question. 

COMPARISON  OF  PREDICTIONS  WITH  OBSERVATIONS 

The  predictions  of  the  vertical  noise  fields  on  the  ground  obtained  from  COMPWR  and 
NOISLAN  can  be  readily  compared  with  data.  From  1958  through  1967,  at  first  the  National  Bureau  of 
Standards  (NBS)  and  after  a  reorganization  .he  Environmental  Science  Services  Administration  (ESSA) 
supervised  a  worldwide  network  of  atmospheric-radio-noise  recording  stations  [30].  Measured  long¬ 
term  noise  data  at  13  kHz  from  12  of  these  stations  have  been  compared  with  COMPWR/NOISLAN 
predictions.  Figures  1  and  2  show  a  comparison  between  the  NBS  experimental  data  and  the  predic¬ 
tions  using  option  II  of  COMPWR  and  options  1  and  II  of  NOISLAN. 

Whereas  many  data  exist  for  vertically  polarized  atmospheric  electric  noise,  few  data  are  available 
for  vertically  polarized  magnetic  noise.  The  Mitre  Corporation  has  made  a  few  nearly  simultaneous 
measurements  of  the  vertically  and  horizontally  polarized  magnetic  noise  fields  at  high  altitudes  using 
loop  antennas.  The  Mitre  Corporation  personnel  used  a  horizontal  loop  antenna  extended  behind  a  jet 
aircraft  in  measuring  the  vertical  magnetic  field  strength.  The  receiver  used  for  the  measurements  was 
an  RCA  VLF  amplitude-probability-distribution  receiver  loaned  by  NRL  for  the  experiments.  A  pair  of 
orthogonal-crossed-loop  antennas  were  used  in  measuring  the  horizontal  magnetic  field.  The  experi¬ 
menters  state  (Ci.  Hirst,  private  communication)  that  their  crossed-loop  system  did  not  have  an  omni¬ 
directional  reception  pattern,  because  of  induced  coupling  from  the  aircraft  structure.  The  crossed-loop 
antennas  were  coupled  to  produce  a  response  to  elliptically  polarized  horizontal  magnetic  noise  fields. 
Neither  the  directions  of  the  major  and  minor  axes  of  the  ellipse  nor  the  direction  of  rotation  of  the 
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Fig.  I  —  Histograms  of  errors  (observed  values 
minus  predicted  values)  in  predicting  vertically  polar¬ 
ized  electric  noise  at  1.1  kHz  using  NOISI.AN  option  I 
and  COMI’WR  option  II 


ERROR  (dB) 


Fig.  2  —  Histogram  of  errors  (observed  values 
'minus  predicted  values)  in  predicting  vertically  po¬ 
larized  electric  noise  at  13  kHz  using  NOISLAN  op¬ 
tion  II  and  COMPWR  option  II. 
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elliptically  polarized  waves  are  available.  These  parameters  are  important  in  comparing  the  data  with 
detailed  theoretical  predictions  from  the  atmospheric-noise  program.  They  were  not  important  to  the 
Mitre  measurement  program.  In  spite  of  these  difficulties,  the  Mitre  data  are  valuable  at  piesent, 
because  they  are  practically  all  the  data  there  are.  The  voltage  outputs  of  the  crossed-loop  antennas 
were  added  together  with  a  90°  phase  shift,  as  is  normally  done  in  a  crossed-loop  system.  The  resultant 
fieu1  strengths  were  then  normalized  to  the  value  of  the  vertical  electric  field  that  an  equivalent  TM 
propagation  wave  would  possess  to  give  the  observed  magnetic-field  value.  These  are  the  values  of  E. 
presented  in  Figs.  3  through  8  and  in  Table  1.  The  values  of^,,  where  %CZ  =  t ?//.,  with  tj  being  the 
impedance  of  free  space,  were  those  measured  on  the  horizontal  loop  antenna. 

From  Figs.  3  through  8  and  from  Table  1,  one  can  note  that  the  equivalent  vertical  electric  noise 
field  Ez  is  predicted  to  be  relatively  constant  from  day  to  night  and  that  it  exceeds  the  equivalent  verti¬ 
cal  magnetic  noise  field  by  about  30  dB  during  the  day  and  by  about  15  dB  during  the  night.  The 
smaller  predicted  value  of  3V.  during  the  day  is  directly  related  to  the  smaller  value  of  the  ionospheric 
reflection  coefficient  h#?j_  that  describes  the  change  of  wave  polarization  from  vertical  electric  to  vertical 
magnetic  (Eq.  (16')). 

An  interesting  case  of  disagreement  with  the  model  is  observed  in  Fig,  4,  where  the  noise  data, 
taken  immediately  after  sunrise  at  the  receiver,  have  a  separation  of  W.  from  E:  more  characteristic  of 
nighttime  than  daytime  conditions.  Such  a  case  could  occur  if  about  50%  of  the  noise  sources  were  in 


10 


"■■'"flftriiespifl 


Location 

. /;.  ~ 

time 

//,  RN 

4S-Noisc  Separation  (dB) 

-  - 

Day 

Nighttime 

17  kHz 

29  kHz 

17  kl 

z 

29  kl 

z 

Opt.  1 

Opt.ll 

Obs. 

Opt.  1 

Opt.ll 

Obs. 

Opt.  1 

Opt.ll 

Obs. 

Opt.  I 

Opt.ll 

Obs. 

47°N,  55°W 

34.5 

35.6 

Eua 

29.5 

28.4 

14.0 

13.5 

— 

12.5 

6.5 

— 

44°N,  65°W 

29.0 

34.5 

— 

26,5 

29.5 

— 

15.0 

15.0 

23.2* 

13.5 

9.5 

13.0* 

38°N,  82°W 

26.0 

32.0 

— 

26.5 

30.5 

— 

12.0 

12.0 

18.1 

13.5 

16.1 

25.0 

31.0 

28.2 

26.0 

29.5 

26.6 

12.0 

13.0 

— 

10.0 

— 

28°N,  80°W 

26.5 

32.0 

35.4 

26.5 

30.0 

26.2 

-2.5 

-3.0 

— 

4.0 

— 

38°N,  84°W 

26.0 

32.0 

— 

26.0 

30.0 

— 

8.0 

8.0 

16.9 

8.5 

7.0 

10.2 

‘Measurements  taken  during  the  transition  from  night  to  day. 


the  nighttime  portion  of  the  globe  and  if  the  effect  of  passage  of  the  nighttime  quasi-TM  waves  through 
the  terminator  region  did  not  suppress  the  fields  but  simply  propagated  them  forward  in  quasi-TE 
modes.  A  future  improved  version  of  the  noise  prediction  model  may  be  constructed  to  take  this 
quasi-TE  conversion  into  account. 

In  general  the  measured  data  are  several  dB  higher  than  the  predictions.  However,  almost  all  the 
E-  measurements  are  within  1  standard  deviation  of  the  predicted  levels  based  on  the  values  of  D„ 
given  in  CCIR  Report  322  (1).  Table  l  presents  a  summary  of  predicted  and  observed  E.  -  $t.  noise 
separations.  The  use  of  option  II  of  the  NOISLAN  program  with  option  II  of  COMPWR  yields  better 
predictions  of  E.  -Sf.  noise  separation  for  daytime  at  =17  kHz  than  does  the  use  of  option  I  of 
NOISLAN  with  option  II  of  COMPWR.  Other  than  for  daytime  at  =17  kHz  the  predictions  are  not 
better  from  use  of  option  II  of  NOISLAN  rather  than  of  option  1.  For  both  options  the  predictions  for 
daytime  fit  the  data  quite  well  but  the  predictions  for  nighttime  do  not  fit  so  well.  The  results  are 
encouraging,  but  more  noise  data  are  needed  to  fully  determine  the  accuracy  of  the  model. 


ADDITIONAL  DEVELOPMENT  OF  THE  MODEL 

Several  possibilities  exist  for  further  development  of  the  atmospheric-noise  model.  Incorporation 
of  TE-mode  propagation  with  conversion  at  the  terminator  may  be  a  helpful  addition.  Also,  the  exten¬ 
sion  of  the  model's  frequency  range  up  to  60  kHz  would  be  useful,  since  the  range  from  30  to  60  kHz 
not  currently  covered  is  of  operational  concern.  Such  an  extension  would  first  require  a  modification  of 
the  COMPWR  and  NOISLAN  programs  to  predict  TM  noise  up  to  60  kHz.  The  propagation  parame¬ 
ters  could  be  found  using  the  WAVEGUID  program  with  an  appropriate  ionospheric  model,  just  the 
way  option  II  of  NOISLAN  was  prepared,  The  validity  of  the  extended  NOISLAN  model  could  be 
checked  against  51 -kHz  vertical  noise  data  collected  by  NBS  [30],  In  preparation  for  this  comparison 
effort,  the  available  51 -kHz  NBS  noise  data  have  been  digitized.  Results  of  the  analysis  may  be  pub¬ 
lished  in  the  future. 
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